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1. Introduction 

There is recently a growing theoretical interest in studying exclusive processes. Exclusive 
production of the Higgs boson is a flag process of special interest and importance. Only a few 
processes have been measured so far at the Tevatron (see pi] and references therein). Khoze, Martin 
and Ryskin developed an approach in the language of off-diagonal unintegrated gluon distributions. 
This approach was applied to exclusive production of Higgs boson [§]. In our recent paper we 
applied the formalism to exclusive production of cc quarks. Quite large cross sections have been 
found ||]. 

The cross section for the Standard Model Higgs production is of the order of 1 fb for Mh = 
120 GeV [^. The dominant bb decay channel is therefore preferential from the point of view of 
statistics. The bb exclusive production was estimated only at higher order |Q]. It was argued that 
the leading-order contribution is rather small using a so-called /j = rule [Q]. Here we show a 
quantitative calculation which goes beyond this simple rule. In our calculation we include exact 
matrix element for massive quarks and the 2 — ;■ 4 phase space. This fully four-body calculation al- 
lows to impose cuts on any kinematical variable one wish to select. Different types of backgrounds 
to Higgs production were studied before e.g. in Ref.[||]. 

2. Formalism 

2.1 The amplitude for pp — ;• ppQQ 

Let us concentrate on the simplest case of the production of qq pair in the color singlet state. 
Color octet state would demand an emission of an extra gluon |Q] which considerably complicates 
the calculations. We do not consider the qqg contribution as it is higher order compared to the one 
considered here. 

In analogy to the Khoze-Martin-Ryskin approach (KMR) [§] for Higgs boson production, we 
write the amphtude of the exclusive diffractive qq pair production pp — ^ p{qq)p in the color singlet 
state as 

''"''HPW2MM) = s.n'^-^^3 J d'qo.Vll'^iquqiMM) 
fgj W,x^ , (7q J ,qif,ti) f°2 [X2,X2,qQ, , ^2,; ) h) 
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%,t1\,tl2,t 



(2.1) 



where A^, Xq are helicities of heavy q and q, respectively. Above /°^ and /2^ are the off-diagonal 
unintegrated gluon distributions in nucleon 1 and 2, respectively. 

The longitudinal momentum fractions of active gluons are calculated based on kinematical 
variables of outgoing quark and antiquark: x\ = ^ exp(+j3) + ^ exp(+y4) and X2 = ^ exp(— 3^3) + 
^exp(— 3^4), where /msj and m^^t are transverse masses of the quark and antiquark, respectively, 
and j3 and y^ are corresponding rapidities. 

The bare amplitude above is subjected to absorption corrections. The absorption corrections 
are taken here in a simple multiplicative form. 
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2.2 gg -^ QQ vertex 

Let us consider the subprocess amplitude for the qq pair production via off-shell gluon-gluon 



fusion. The vertex factor V^ '^^ = V^'^ {qi,q2,ki,k2)in expression ( 2. 1 ) is the production amplitude 
of a pair of massive quark q and antiquark q with helicities A^, Xq and momenta k\,k2, respectively. 
The color singlet qq pair production amplitude can be written as 



\aJ('?1'^2,^i,^2) = «^«v\'a^ {quqiMM), (2.2) 

The tensorial part of the amplitude reads: 

T7/JV / 1 1 \ 2- n \( V q\-k\-m „ y qi-k2+m ^\ 

(2.3) 

The coupling constants g1 — )• gsil^-^i )8s{l^r2)- ^^ ^^ present calculation we take the renormalization 
scale to be n^^ = IJ.^2 = ^gs/^ or ^qs- The exact matrix element is calculated numerically. 

2.3 Off-diagonal unintegrated gluon distributions 

In the KMR approach the off-diagonal parton distributions (i=l,2) are calculated as 

.KMR. ^2 ,,2 ,x „ d[g{xi,kf)Si/2{kf,^^)] 



R, 



dg{Xi,kfj^ . c /^2 ,.2 



dlogkf 



'-2^\k2=QlS,;2{Qt„^')F{ti), (2.4) 



where 51/2 (^?iM^) is a Sudakov-like form factor relevant for the case under consideration. The 
last approximate equalities come from the fact that in the region under consideration the Sudakov- 
like form factors are somewhat slower functions of transverse momenta than the coUinear gluon 
distributions. While reasonable for an estimate of gluon distribution it may be not sufficient for 
precise calculation of the cross section. It is reasonable to take a running (factorization) scale as: 
W=M|=^y4orM2^-. 

The factor Rg here cannot be calculated from first principles in the most general case of off- 
diagonal UGDFs. It can be estimated in the case of off-diagonal collinear PDFs when x' <^x and 
xg = x^^{l —x)". Then Rg = ^7f~ nA+4) Typically Rg ^ 1.3 - 1.4 at the Tevatron energy. The 
off-diagonal form factors are parametrized here as F(f ) = exp (Sofff ). In practical calculations we 
take Boff = 2 GeV^^. In the original KMR approach the following prescription for the effective 
transverse momentum is taken: Q^ , = min ( '/o p ^i r ) ^^^ Q\ t ~ ^^ ( ^o r ' ^2 ? ) • ^^ evaluating /i 
and /2 needed for calculating the amplitude (Ol) we use different collinear distributions. It was 



proposed |g|| to express the 5 1/2 form factors in Eq. ( |2.4[ ) through the standard Sudakov form factors 
as: 

,2 ,,2\ 



S,/2{qf,^') = jTg{qj,^^). (2.5) 
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2.4 Cross section 

The cross section is calculated as 

da = —\.y^2^A\^{2n,f5'^{pa+Pb-p\-p2-p-i-pA 



(Ppi 



d^P2 d^p3 d^p4 



{InyiEi {271Y2E2 {InyiEs {InyiE^ ' 

(2.6) 

The details how to conveniently reduce the number of kinematical integration variables are given 

elsewhere. 



3. Results 



3.1 /?p — )• ppcc 

Let us proceed now with the presentation of differential distributions of charm quarks produced 
in the EDD mechanism. In this case we have fixed the scale of the Sudakov form factor to be 
/I = Mf-cjl. Such a choice of the scale leads to a strong damping of the cases with large rapidity 
gaps between q and q. 

In the left panel of Fig. [T] we show distribution in rapidity. The results obtained with the KMR 
method are shown together with inclusive gluon-gluon contribution. The effect of absorption leads 
to a damping of the cross section by an energy-dependent factor. For the Tevatron this factor is 
about 0.1. If the extra factor is taken into account the EDD contribution is of the order of 1 % of the 
dominant gluon-gluon fusion contribution. 

The corresponding rapidity-integrated cross section at ^fs = 1960 GeV is: 6.6 jxh for the exact 



formula, 2.4 /ib for the simplified formula (see Eq. (|2.4[)). For comparison the inclusive cross 
section (gluon-gluon component only) is 807 jxh. 
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Figure 1: Rapidity distribution of c or c (left) and transverse momentum distribution of c or c (right). 
The top curve is for inclusive production in the A:, -factorization approach with the Kwiecinski UGDF and 
ji^ = 4m^, while the two lower lines are for the EDD mechanism for the KMR UGDF with leading-order 



collinear gluon distribution [|11|. The solid line is calculated from the exact formula and the dashed line for 
the simplified formula (when only derivative of collinear GDF is taken). An extra cut on the momenta in the 
loop QJcut = 0.26 GeV^ was imposed. Absorption effects were included approximately by multiplying the 
cross section by the gap survival factor Sg = 0.1. 
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In the right panel of Fig. |I| we show the differential cross section in transverse momentum 
of the charm quark. Compared to the inclusive case, the exclusive contribution falls significantly 
faster with transverse momentum than in the inclusive case. 

In Fig. ^ we show the distribution in the invariant mass of the c c system. Compared to the 
inclusive case the invariant mass distribution for the FDD component is significantly steeper. This 
is due to the Sudakov-like form factor which, according to the procedure described above, damps 
the cross section for large invariant masses. 
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Figure 2: Invariant mass distribution of the cc pair (left) and distribution in the transverse momentum of the 
cc pair (right). The other details are the same as in the previous figure 



As in the inclusive case within the ^, -factorisation approach the cc pair possesses the transverse 
momentum different from zero. The corresponding distribution is shown in the right panel of 
Fig. 1^ The distribution for the exclusive case (the two lower lines) is much narrower compared to 
the inclusive case (the upper line). 

3.2 PP — 7- ppbb 

In parallel to the exclusive bb production, we calculate the differential cross sections for ex- 
clusive Higgs boson production. Compared to the standard KMR approach here we calculate the 
amplitude with the hard subprocess g*g* — )• H taking into account off-shellness of the active glu- 
ons. The details of the off-shell matrix element can be found in Ref. ||^. In contrast to the exclusive 
Xc production [^, due to a large factorization scale ~ Mh the off-shell effects for g*g* — ;■ H give 
only a few percents to the final result. 

The same unintegrated gluon distributions based on the coUinear distributions are used for the 
Higgs and continuum bb production. In the case of exclusive Higgs production we calculate the 
four-dimensional distribution in the standard kinematical variables: y,ti,t2 and <p. Assuming the 
full coverage for outgoing protons we construct the two-dimensional distributions do /dycf'pt in 
Higgs rapidity and transverse momentum. The distribution is used then in a simple Monte Carlo 
code which includes the Higgs boson decay into the bb channel. It is checked subsequently whether 
b and b enter into the pseudorapidity region spanned by the central detector. 

In Fig. H we show the most essential distribution in the invariant mass of the centrally produced 
bb pair, which is also being the missing mass of the two outgoing protons. In this calculation we 
have taken into account typical detector limitations in rapidity —2.5 <yh,yi<2.5. We show results 
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Figure 3: The bb invariant mass distribution for y/s = 14 TeV and for —2.5 < y < 2.5 corresponding to the 
ATLAS/CMS detectors. The absorption effects were taken into account by muhiplying by the gap survival 
factor 5g = 0.03. 
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Figure 4: The bb invariant mass distribution for y/s = 14 TeV and for b and b jets in the rapidity interval 
—2.5 <y < 2.5 corresponding to the ATLAS detector The absorption effects for the Higgs boson and the 
background were taken into account by multiplying by the gap survival factor Sq = 0.03. The left panel 
shows purely theoretical predictions, while the right panel includes experimental effects due to experimental 
uncertainty in invariant mass measurement. 



with different coUinear gluon distributions from the literature: GRV [|11[], CTEQ []12|], GJR [ |13| ] 
and MSTW [|l^]. The results obtained with radiatively generated gluon distributions (GRV, GJR) 
allow to use low values of Qt = qot,qit,q2t whereas for other gluon distributions an upper cut on Qt 
is necessary. The integrated double-diffractive bb contribution calculated here seems bigger than 
the contribution of the exclusive photoproduction of bb estimated in [15|. The lowest curve in Fig.|| 
represents the 77 contribution [J8|]. While the integrated over phase space 77 contribution is rather 
small, it is significant compared to the double-diffractive component at large M^g > 100 GeV. This 
can be understood by a damping of the double diffractive component at large M^^ by the Sudakov 
form factor [Q, |3|]. In addition, in contrast to the double-diffractive component the absorption for 
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Figure 5: The bb invariant mass distribution for y/s = 14 TeV for a limited range of b and b rapidities: 

-l<y<l. 



the YY component is very small and in practice can be neglected. 

In the left panel of Fig.0 we show the double diffractive contribution for a selected (CTEQ6 



[]12p) collinear gluon distribution and the contribution from the decay of the Higgs boson including 
natural decay width calculated as in Ref. [|^], see the sharp peak at M^^ = 120 GeV. The phase 
space integrated cross section for the Higgs production, including absorption effects with So = 0.03 
is somewhat less than 1 fb. The result shown in Fig.^ includes also the branching fraction for 
BR(// — )• bb) « 0.8 and the rapidity restrictions. The second much broader Breit-Wigner type 
peak corresponds to the exclusive production of the Z° boson with the cross section calculated as 
in Ref. [p^]. The exclusive cross section for ^/s = 14 TeV is 16.61 fb including absorption. The 
branching fraction BR(Z*' — )• bb) ss 0.15 has been included in addition. In contrast to the Higgs case 
the absorption effects for the Z^ production are much smaller [|T7[]. The sharp peak corresponding 
to the Higgs boson clearly sticks above the background. In the above calculations we have assumed 
an ideal no-error measurement. 

In reality the situation is, however, much worse as both protons and in particular b and ^ jets are 
measured with a certain precision which automatically leads to a smearing in M^g . Experimentally 
instead of M^g one will measure rather two-proton missing mass (Mpp). The experimental effects 
are included in the simplest way by a convolution of the theoretical distributions with the Gaussian 
smearing function G{M) = ^i-^exp ( ^ '^^"' ] with a = 2 GeV [[T^, |l9|] which is determined 



iTia 



mainly by the precision of measuring forward protons. In the right panel we show the two-proton 
missinging mass distribution when the smearing is included. Now the bump corresponding to the 
Higgs boson is below the bb background. With the experimental resolution assumed above the 
identification of the Standard Model Higgs seems rather difficult. The situation for some scenarios 
beyond the Standard Model may be better [^, 21]. 

Can the situation be improved by imposing further cuts? In Fig. || (left panel) we show the 
result for a more limited range of b and b rapidity, i.e. not making use of the whole coverage of 
the main LHC detectors. Here we omit the Z" contribution and concentrate solely on the Higgs 
signal. Now the signal-to-background ratio is somewhat improved. This would be obviously at the 



Exclusive production of heavy quark pairs Antoni Szczurek 



expense of a deteriorated statistics. Similar improvements of the signal-to-background ratio can be 
obtained by imposing cuts on jet transverse momenta. Detailed studies of the role of cuts will be 
discussed in [^. 

We are indebted to Valery Khoze, Misha Ryskin, Andy Pilkington and Christophe Royon for 
a discussion and exchange of useful information. 
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